Numerical Renormalization Group at Criticality 
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Abstract 

We apply a recently developed numerical renormalization group, the corner- 
transfer-matrix renormalization group (CTMRG), to 2D classical lattice models 
at their critical temperatures. It is shown that the combination of CTMRG and 
the finite-size scaling analysis gives two independent critical exponents. 
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The renormalization group is a basic concept in statistical physics ||T|, |2|. The 
real-space renormalization group (RSRG) has been applied widely to critical phe- 
nomena fl], g. A resent progress in RSRG is the development of the density matrix 
renormalization group (DMRG) by White |, |]. The DMRG was originally formu- 
lated for one-dimensional (ID) quantum models, and was shown to be applicable 
to 2D classical systems as well|^. Recently, Ostlund and Rommer found a vari- 
ational principle hiding behind DMRG 0]. They showed that the ground state 
wave function obtained by DMRG is in a form of a matrix product ^, This 
variational background enables us to accelerate the numerical calculation of DMRG 
Quite recently Martin-Delgado and Sierra have obtained a unified analytic 



formulation of conventional RSRG and DMRG [^, 13] 



Baxter established another RSRG for 2D classical systems by using the corner 



transfer matrix (CTM) |14, 15, 16]. His method is known as a natural extension 



of variational methods, such as Kramers- Wannier approximation |17], Kikuchi's ap- 
proximation ]|l^ , and the cluster variational method ]19]. It should be noted that 
Baxter's method and White's DMRG method have many features in common. In 
particular, both are RSRG based on a variational principle. On the basis of this 
fact, Nishino and Okunishi have formulated a new renormalization group procedure 
for 2D classical lattice models |2^, which we will refer 'corner-transfer-matrix 
renormalization group (CTMRG)' in the following. It has been confirmed that the 
CTMRG precisely determines thermodynamic functions of 2D classical systems in 
the thermodynamic limit. For example, calculated internal energy of the square- 
lattice Ising model is -0.70704 at Tc [20], whose numerical precision is comparable 
to a resent Monte Cairo result ]21]. 

In this letter, we show that the CTMRG method determines not only the one- 
point functions, but also critical exponents of 2D classical spin systems. We apply 
CTMRG to finite size systems at their critical temperatures, and calculate one-point 
functions at the center of the system. The finite size corrections to these one-point 
functions are found to obey the finite-size scaling (FSS) behaviors determined by 



two independent critical exponents ]^, |23[ ]. 

The CTMRG is a numerical method which can evaluate the partition function of 
finite size systems. Consider a square cluster of a classical spin system whose linear 
dimension N is an odd integer. The cluster consists of four sub-clusters of the size 
{N + l)/2, which are named 'corners' according to Baxter; The partition function 

is calculated as the trace of a matrix P[n +i)/2 = (^(Af+i)/2) i where ^(Ar+i)/2 is so 
called 'corner transfer matrix' (CTM) which transfers column spins into row spins 
]16]. It, however, is difficult to deal with ^(Ar+i)/2 exactly for large N because 
the dimension of the matrix increases rapidly with N . The point of CTMRG is, 
that the matrix p is identical to the density matrix that appears in White's DMRG 
]^, ^. This correspondence enables us to transform A(^^^iy2 into a renormalized 
one ^(Af+i)/2- We determine the dimension m of the reduced matrix ^(Ar+i)/2 so 
that we can deal with 74(jv+i)/2 realistic numerical calculations. By using a 
recursive relation between ^(Ar+i)/2 and ^(Ar+3)/2; we can increase the linear size of 
the corner one by one, and obtain ^(Ar+i)/2 from = 3 to oo one after another; this 
recursive procedure corresponds to the 'infinite chain method' in DMRG ]Q, ^. The 
renormalized CTM ^(Ar+i)/2 thus obtained gives the approximate partition function 
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as Zn = Tr^(jv+i)/2- Because of the variational nature of the CTMRG, Zn gives a 
lower-bound for the exact partition function Z^r = TrA(7v+i)/2- 

The error in the partition function dZ^ = Zj^j—Z^ is related to two characteristic 
length scales. In what follows, we denote the largest and the second largest eigenval- 
ues of the row-to-row transfer matrix as Ao(A^, m) and Ai(A^, m) respectively, which 
can be calculated by CTMRG. The first characteristic scale is the correct correlation 
length for a finite system defined as follows: 

Ao(iV,oo) 

According to the finite-size scaling theory, is a decreasing function of A'" as 

~ at the critical temperature T = Tc- The other scale is an effective correlation 
length ^{m) for finite m in the thermodynamic limit: 

1/^ m = log — — -, (2) 

Ai(c>o, m) 

This latter scale appears due to the restriction imposed on the size of the matrix 
^(Ar+i)/2- As long as the condition <C iim) is satisfied, Z;^ is a good approxima- 
tion for Zat. But in the opposite case ^ ^ijn)-, will be appreciably smaller 
than Z]\f. Therefore, we expect to observe crossover from the correct system-size 
dependence of the partition function to finite-m behavior. For off-critical case, we 
can make the condition <C ^(m) satisfied by taking large enough m |2^. For the 
critical case, on the other hand, we have to deal with the above crossover properly. 

In addition to the partition function, we can also calculate the local energy 
E{N) and the order parameter M{N) at the center of the square clusters of the 



linear size A^[20|. According to the above discussion on the crossover, we expect 
that the following two-parameter finite-size scaling form |2^, ^] holds for these 
thermodynamic functions at T^- The N dependence of E{N) is 

E{N) - Ei^) = N'/'^-^fiam)/N), (3) 

where is one of the critical exponents (the correlation length exponent) and d = 2 
is the spatial dimensionality. In addition to the leading finite-size scaling form ~ 
j\ji/u-d^ we introduced an unknown scaling function / which describes the crossover. 
The asymptotic form of the scaling function will be f{x — > 0) ~ rj.-'^/'^+d g^j^j _^ 
oo) ~ const.. Similarly, the size dependence of the local order parameter M{N) is 
expected to be 

M{N) = N-'^'^-^+'i'>/^g{C{m)/N) (4) 

with another set of critical exponent rj (the anomalous dimension of the spin) and 
a scaling function g introduced; The asymptotic behavior of g will be g{x — > 0) ~ 
^{d-2+r})/2 g^^^ _^ ^ const.. Assuming these two scaling form, we can 
estimate two independent critical exponents i' and rj as well as E{oo) from the data 
given by the CTMRG. 

As examples of the 2D classical lattice models, we deal with Ising model and 
3-state Potts model on the square lattice. The exact values of the critical exponents 
appearing in the scaling forms are v = \ and = 1/4 for the Ising model and 
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u = 6/5 and r] = 4/15 for the Potts model. In order to calculate M{N), we impose 
the ferromagnetic boundary conditions, that is, all the spins at the boundary of the 
clusters take the same value. 

Figure 1 shows the local energy E{N) = {aiai+i) of the Ising model against 
1/N at Tc = 2.269185314. We set m = 148; the result for m = 4 is also shown 
for comparison. The linear dependence of E{N) on 1/A^ is observed in a wide 
range of 1/N, which is consistent with the leading FSS behavior in Eq. (3) with 
u = 1. In fact, the least-square fitting of the data in the range 21 < < 401 
to Eq. (3) gives u = 0.993 and E{oo) = 0.70704, which is close to the exact value 
£^(oo) = 1/V2 = 0.70711. Figure 2 shows the spin polarization M{N) = {a) for 
m = 4, 10, and 148. The N dependence is well expressed by N~^/^ as we expected. 
The crossover effect discussed above is clearly seen in the figure: M{N) for different 
m deviate from N~^/^ behavior one by one when N increases. By the least-square 
fitting in the range 21 < A'' < 401 wc get r] = 0.2506. 

Now we check the full two-parameter scaling form for E{N) taking the crossover 
effect expressed by the scaling function / into account. In figure 3 we plot N {E{N) — E{oo)} 
against N/^{m) for m = 4 ~ 13; we choose relatively small m in order to observe 
small ^{m)/N region. All the data are really collapse into a single curve, and thus 
the two-parameter scaling assumption, Eq. (3), is justified. 

Figures 4 and 5 show the local energy E{N) = {5{aiai^i)) and the local or- 
der parameter M{N) = (5(0, cr)), respectively, of the 3-state Potts model at Tc = 
1.989944809. Again the expected A^ dependence, E{N) - E{oo) ~ A^-^/s and 
M{N) ~ N'"^/^^, are observed. The calculated critical exponent u is 0.830 in the 
range 5 < A?^ < 31, and is 0.809 in the range 4 < A?^ < 201. For the exponent 77 we 
get 0.266 in the range 5 < A/" < 31, and 0.267 in the range 5<N< 201. 

The obtained exponents u and 77 of the 3-state Potts model seem to have lower 
accuracy than those of the Ising model. The reason may be attributed to that we 
used the same value of m for both models. Since the Potts model has a larger spin 
degree of freedom than the Ising model, a larger size of matrix may be required to 
retain the same order of accuracy as the Ising case. 

The authors would like to express their sincere thanks to Y. Akutsu for valuable 
discussions. T. N. thank to S. R. White for helpful discussions about the DMRG 
method. He is also grateful to G. Sierra for discussions about RSRG. The present 
work is partially supported by a Grant-in- Aid from Ministry of Education, Science 
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Figure Captions 

Fig. 1. Nearest neighbor spin correlation function of the Ising model. 

Fig. 2. Spin polarization of the Ising model. 

Fig. 3. Two-parameter scaling function /(^(m)/iV). 

Fig. 4. Nearest neighbor spin correlation function of the 3-state Potts model. 

Fig. 5. Spin polarization of the 3-statc Potts model. 
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